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Abstract 



Our main goal in this paper is to extend to any system of coupled quadratic Hamiltonians 
some properties known for systems of quantum harmonic oscillators related with the Brownian 
Quantum Motion model. In a first part we get a rather general formula for the purity (or the 
, linear entropy) in a short time approximation. In a second part we establish a master equation 

qq ■ (or a Fokker-Planck type equation) for the time evolution of the reduced matrix density for 

bilinearly coupled quadratic Hamiltonians. The Hamiltonians and the bilinear coupling can be 
time dependent. 

Moreover we give an explicit formula for the solution of this master equation so that the time 
evolution of the reduced density at time t is connected with the reduced density at initial time to 
for to < t < to + t c where t c G]0, do] is a critical time but reversibility is lost for t > to + t c . 

-I— > 

1 Introduction 

The general setting considered here is a quantum system (S) interacting with an environment (E) . The 
' total system (S)U(E) is supposed to be an isolated quantum system and we are interested in dynamical 

properties of (S) alone, which is an open system because of its interactions with (E). In particular 
during the time evolution the energy of (S) is not preserved and its evolution is not determined by a 
Schrodinger or Liouville-von Neumann equation unlike for the total system (S) U (E). 
The Hilbert space of the total system (S)U(E) is the tensor product % = 'Hs^'H.eQ and its Hamiltonian 
H is decomposed as follows: 

H = H S ®I + I®H E +H I = H + V (1.1) 

s v ' 

Ha 

(S)U(E) is isolated by assumption and its evolution obeys the Schrodinger (or Liouville-von Neumann) 
; i " equation with the Hamiltonian H. 

All Hamiltonians here are self -adjoint operators on their natural domain, H and Hj — V are defined 
in W, Hs m Hsi He hi He- Moreover most of the results stated here are valid when the Hamiltonians 
are time dependent, assuming that their propagators exist as unitary operators in the corresponding 
Hilbert space. 

Quantum observables are denoted with a hat accent, the corresponding classical observables (also 
named Wigner function or Weyl-symbols) are written by erasing the hat. 

We don't give here more details concerning the domains, these will be clear in the applications. 
We assume that the interacting potential V has the following form 

V= ]T SjQEj (1.2) 

l<j<m 



where Sj and Ej are self-adjoint operators in Tisj He respectively. 



1 when the spaces are infinite dimensional we mean here that Hs 8Mb is the Hilbert tensor product, completion of 
the algebraic tensor product 
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Recall that a density matrix p is a positive class trace operator with trace one (trp =1). p is a 
state of the total system (5) U (E) in the Hilbert space W. The time evolution of p obeys the following 
Liouville-von Neumann equation: 

p = i-\H,pl P :t^p{t) (1.3) 

Time derivatives are denoted by a dot, [•, •] denotes the commutator of two observables. 
We assume that for t = the system and the environment are decoupled : 

p(0) = ps(0) <g> Ae(0), ps,e(0) are density matrices in Hs,e (1.4) 

and the state /Ss(0) of the system is pure i.e is an orthogonal projector on a unit vector tp of Hs- A 
non pure state will be called a mixed state. A density matrix ps is a mixed state if and only if ps has 
an eigenvalue A, such that < A < 1. 

Si, if p is a pure state then p = ±1^, tp e ||^| = 1 where n^(r/) = (tjj,r])ip. 

The general problem for open systems is to describe the evolution ps(t) of the density matrix of the 
system (S). In particular an important physical question is: when the state ps(t) is pure or mixed? 
To decide if a state is pure or not we consider the purity function 

Pur(t) = tr s (ps(t) 2 ) 

The function Sg(t) — 1 — p ur (t) is called the linear entropy. The Boltzman-von Neumann entropy is 
Sbn = -trs(ps^ogps)- 

We clearly have < p U r{t) < 1. It is not difficult to see that ps(t) is pure if and only if p ur (t) = 1. 
With our assumption, is a maximum for p ur so we have iw(0) < and if p ur (0) < then for 
< t < e, ps(t) is a mixed state for e small enough. 

For isolated (closed) systems, if the initial state is pure then it stays pure at every time, as it can be 
easily seen using equation (|1.3[) . 

So our first step in this paper is to compute p ur (0). In particular for a large class of models, including 
the Quantum Brownian Motion, we shall prove that p ur (0) < so for these models the state of the 
system becomes very quickly decoherent and entangled with the environment. 

Our second goal is to compute the time evolution ps{t) of the state of the system. We know that 
ps (t) does not satisfy a Liouville- von Neumann equation because the system is open. 
We shall prove that for any time dependent quadratic systems, ps{t) obeys an exact "master equation" 
similar to a Fokker-Planck type equation, with time dependent coefficients, and from this equation 
we can get an explicit formula for ps(t)- 

This could be a first step for more general systems in the semi-classical regime . 

Let us recall a mathematical definition for ps(t) and related properties and notations. 

This can be done by introducing partial trace (or relative trace) over the environment for a state p of 

the global system. 

Definition 1.1 Let A be a class-trace operator in Us ®T-Le- We denote tr^A the unique trace-class 
operator on %s satisfying, for every bounded operator B on Us, 

tr s ((tr B i)s) - tr (A[B ® I E )) (1.5) 

We have denoted "tr" the trace in the total space H, trs,E the trace in %s,e respectively. 

Of course we have ti{B <g> C) = B ■ (tr^C) if B and C are trace class operators in Hs and He 

respectively. 

If p is a density matrix in T-L then tr^p is a density matrix in Hs, called the reduced density matrix 
or the reduced state. 

From (|1.5[) we can compute the matrix element of trEA by the formula 

(tp, {tr E A)<p) = tr (A(II^, v <gs I E j) , W>, <p e Hs- 
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The partial trace was introduced in quantum mechanics to explain quantum phenomena like entan- 
glement and decoherence (see [TSJ [T^] for more details). 

Let G Hs®He, ^ = J]] ipj <& r)j. $ is a pure state in H = Hs %He- Let us compute the partial 

l<j<JV 

trace of 11$ in He- Applying Definition ll.il we easily get 

tr B (n*)= (VjiVk)^,^ (1.6) 

l<j,fc<JV 

where IL^.^ is the rank one operator in Hs ■ ^i^.^ k {(p) = (^pj,(p)tpk- 

Decoherence means in particular that there exists a basis {rjj} such that K%)%)| becomes very small 
for j 7^ k, so that trg(n$) is very close to H^ hence quantum interferences for the system (S) 

i<j<N 

are lost. 

Formula (|1.5[) is an operator version of Fubini integration theorem as we can see in the Weyl quanti- 
zation setting. 

Let H s = L 2 (R d ) , He = L 2 (R N ). Denote z = (x,f) G R 2d , u = (y,rf) G M. 2N , A = ct^A the Weyl 
symbol of A in the Schwartz space S(HL 2 ( d+N ') (see [6] for more details) then the Weyl symbol of IteA 
is 

a w {tv E A){z) = (2tt)- n [ A{z,u)du. (1.7) 
Our main applications here concerns the Weyl-quadratic case where : 

• Hs, He are quadratic Hamiltonians (generalized harmonic oscillators) respectively in the phase 
spaces R 2d , M. 2N . 

• Hs, He are quantum Hamiltonians (Weyl quantization of Hs respectively He) in L 2 (M. d ) re- 
spectively L 2 (R N ). 

• ps is the projector on tJj G S(R d ). 

• pe is a Gaussian with means. 

• Sj resp. Ej are linear forms on M. 2d resp. K 2Ar (bilinear coupling). 

In applications one consider thermal equilibrium states for environment: pe — Z((3)~ 1 e~^ KE where 
Ke is a quadratic form, positive-definite and j3 = A, T > is the temperature, Z(/3) = tr(e~@ KE ). 
In this setting it is possible to find the exact time dependent master equation satisfied by ps(t), as we 
shall see in the second part of this paper {pe can be any environment state) and to solve this equation 
when the classical dynamics of the total system (S) U (E) is known. 

Let us remark that there exists a general formula, called Kraus formula, giving the time evolution 
of the reduced system: p s {t) = tr E (U(t)p(Q)U* (t)) where p(0) = ps (0) ® p E (0) and U{t) = e- %tH (see 
[HIE]): _ 

ps(t) = J2K j , i (t)p s (0)K j M* (1-8) 
it 

where the Kraus operators Kj^(t) depends on pe{0 and U(t). But the operators Kj^(t) are not 
explicitly known, so formula (|1.8p is not easy to use to get properties of ps(t)- 

2 A general computation 

It is easier here to consider the interaction formulation of quantum mechanics to "eliminate" the "free" 
(non-interacting) evolution: Uo(t) — e~~ ltH ° . In this section we consider the general setting described 
in the introduction. No restrictions are necessary for the Hamiltonians Hs, He, we only need that 
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the interaction Hamiltonian has the shape (jl.2|) . As it is well known the time evolution (|1.3j) is given 
by the Schrodinger evolution 



p(t) = U(t)p(0p(t)*, U(t)=c 



-itH 



For simplicity we assume here that all Hamiltonians are time independent but the results are also 
valid if H, Sj,Ej are time dependent. In this case we have U(t) = U (t, 0) where is the initial time. 
Introduce the inter-acting evolution 

pW(t) = u (typ(t)u (t), 

(|1.3p becomes 

p I) =i- 1 W(t),n (2-9) 

where 

V(t) = U (t)*VU (t) = S 3 (t)(g>E 3 (t), 

l<j<m 

with Sj(t) = e lt " s S.je- ltfls and E 3 {t) = c UHe E 3 e~ ltHE . The useful fact here is that the shape of the 
interaction V(t) is the same for all t G M (|1.2p . Let us remark again that technical assumptions are 
necessary for rigorous statements. 
We see easily that 

p { P(t) := tr E (p^)(t)) = c^p s (t)e~^. 
We want to consider the purity p ur (t) = trsipgit) 2 ). We shall use now the simpler notation f(t) = 

Using cyclicity of the trace we have 

p ur (t) = 2tr s {f{t)r(t)) (2.10) 

p ur (t) = 2(tr s (f(^'(<))+tr s (^)) 2 ) (2.H) 

We have, at time t = 0, 

r(0) = i-hr E [V,p s ®p E ], (2.12) 

Ho) = r l tx E [ir,ps®PE]-txE[V,[V,ps®PE]]. (2.13) 



Lemma 2.1 Assume that ps(0) is a pure state and that tr E (p E (0)E 3 ) = for every 1 < j < ni. 
Then we have 

r (0) = 0. 

Proof. For simplicity we denote p E — p E (0), ps — Psity- Using the splitting assumption on V we 
have 

[V, ps <8> Pe] = [ Y, S j ®E j ,p s ®pE] (2.14) 
Y Sjps^EjpE- psSj^psEj (2.15) 

l<j<m 



So taking the .E-trace we have 



ti E ([V,ps O Pe]) = Y [ s iiPs]toE{PEEj 

1< 7 <m 



□ 
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Lemma 2.2 We have 

tTs(tr E [V,ps®PE]ps)=0 (2.16) 

Proof. This follows from the definition of the relative trace ti E and cyclicity of the trace. 
We have 

trs (tr E [V, p s ® Ae]AsJ = tr f [F, ps ® Pe]ps ® l) 

This is because ps® Pe and ps ® 1 commute. □ 
Finally, using the identity 

[F, [V, p}} = V 2 p + pV 2 - 2VpV 
with the same argument as in Lemma |2. II we get the following proposition: 
Proposition 2.3 Under the condition of Lemma \2.1\ we have 

PuAO) = -4tr (Vp(p)V(p£ <g> I)) (2.17) 

with p(0) = p s (0) ® p B (0) and ^ = I - p s (0). 

In particular if (pg €5 1) VpV" is not in the range of p-g (8 1 t/ien p ur (0) < and /or any t ^ 0, t small 
enough, ps{t) is a mixed state. 

The formula (|2.17[) can be written in a more suggestive form, using the decomposition of V to separate 
the contributions of the system and its environment. The result is 

p ur (0) = -4 ]T (tr s (p S S r S 3 ) - ixstfsSypsSjj) tv E (p E E r E 3 ) (2.18) 

l<j,j'<m 

This formula has the following statistical interpretation, introducing the quantum covariance matrices: 
rjf} = tr s (ps4-'4-) " tes(ps§j>Ps§j), rfj - ME^PEEfEi) (2.19) 

so we have 

p„ r (0) = -4tr Rm (r( s )(r( £ )) T ). (2.20) 
Then we get a short time asymptotic expansion using Taylor formula 

Pur(t) = ^^t 2 +0(t 3 ), t\0. 

So if from formula (I2.20p we can infer that p ur (0) < then we get that ps(t) is a mixed state for t > 
small enough. 

3 Application to the quadratic case 

We consider here the Weyl-quadratic case described in the introduction (all Hamiltonians are now 
supposed to be quadratic). 

We use the following basic property of Weyl quantization in R d : If A, B are Weyl symbols such that 
B is in the Schwartz space iS(M d ) and A is a polynomial (or like a polynomial) then we have 

tr i2(R d)(AB) = (27r)~ d / A(z)B{z)dz. 

Let us compute p ur (0) using formula (|2.19|) and the Weyl symbols. 

Note that 2ir~ d ps(z)dz and (2tt)~~ n p E (u)du are probabilities laws denoted by irs{z)dz and ir E (u)du 
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respectively. 

Ps = /3(0) is a projector on a state ipo of %s so we can compute 
trs(psSfpsSj) = (ipo, Sjtp^iipo, Sj'ipo) 

= (/ ns(z)S j (z)dz)( f 7T S {z)S f {z)dz). (3.21) 



Let us introduce the commutators 

s j'j = [Sj'iSj], efj = [Ej',Ej] 



We have 



tes{psSj>Sj) = / n s (z)S r (z)S j (z)dz + — s fj 



1 

tiEipEEfEj) = I ■K E {u)E jl (u)E j (u)du+^ : e rj . (3.22) 



Then we get 



Pur(Q) = -4V f / nsWSfWSjWdz 

77' UR2d 

- ( / 7r s (z)5 J (z)dz)( / ir s (z)S j >(z)dz) ) ■ V>y/ .,> 

Let us introduce the classical (symmetric) covariance matrices 

Eg? = /" KaWSfWSjWdz - ( / 7rs(z)S J -(«)ds)( / 7rg(z)^( z )dz) 

and 

So we get 



(3.23) 



s iy = / ^^E^Ef^E^du. 



p ur (0) = -4tr K2m (S( s )s( £ ))+ £ S/J e^, (3.24) 

We remark that in formula (|3.24[) r.h.s the first term is classical and the second term is a quantum 
correction q c which vanishes if the coupling between the system and its environment is only in position 
(or momentum) variables ( as it is usually in the literature). 

Assume now that Sj(z) — Zj and Ej(u) = Gj.iui, where 1 < j < 2d. Gj t i is a 2d x 2N 

l<i<2N 

matrix with real coefficients. With these notations we have for the quantum correction 

q c = 2_j Sj'jej'j =tT R 2d(J s GJ E G T ) (3.25) 

l<j,j'<2d 

where Js, Je are respectively the matrix of the canonical symplectic form in M. 2d and R 2JV : crs(z, z') — 
z ■ Jsz', <te{u,u') = u ■ Jeu' (the scalar products are denoted by a -. 
For the "classical part " we have 

- 4tr R 2 m (£( s >£( £; )) = -4tr(G T Cov ps GCov p J (3.26) 

where G T is the transposed of the matrix G. 
So we have proved the formula: 



Pur(0) = -4tr R 2«(G T Cov 7rs GCov TB ) + tY K , d {J s GJ E G T ) (3.27) 
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where Cov^ is the covariance matrix for the quasi-probability tt (the Wigner function of a pure state 
ip is non negative if and only if ip is a Gaussian by the Hudson theorem [5]) 

In the physicist literature [SJ [3J [TU] the following example is often considered: a quantum oscil- 
lator system is coupled with a bath of oscillators in a thermal equilibrium by their position vari- 
ables. This is known as the Quantum Brownian Motion model. So we choose d = 1, ps = Pt,e, 
V(z,u) = 



l<j<N 



ci, • • • , cat are real numbers, z = (x, £) £ Mr, u = (y, 77) G R and pe — Pt,e where 

p TiE = (2t) n exp I t ^ '• 

\ 1<3<2N 

t is related to the temperature T by the formula r = tanh(T _1 ) (taking the Boltzmann constant 
k = 1). So we have < r < 1. 
With these parameters we get 



Pur(0) = (Cj 

r 



+ c 2 N )Cov x/x ('4> ) 



(3.28) 



where Cov x ^ x (ip ) = J R x 2 \tpo(x)\ 2 dx — (J R x\tpo(x)\ 2 dx) 2 . If the oscillators of the bath have different 
temperatures T\ , • • • , TV then the result is 







p ur (0) = -2(^i + • 


• + ^)Cov x . x (Vo) 







where Tj = tanh(T J 1 ). 

A more general bilinear coupling is 



\l<j<Ar l<j<Ar J \l<j<N l<j<N 

The quantum correction q c is 

g c = 2(d 1 c N+1 H d N c 2N ) - 2(cid N+ i H c N d 2N ) 

Assume for simplicity that Cj = for j > N+2 and dj = for j > 2. Then we find after computations: 



Pur{0) = ( c i H h cir + c Ar+1 )var(-0o) 



+2cidiCov a! j(V'o) + c^Cov^^o) + 2cAr + id a 



(3.29) 



So we see that new quantum corrections appear when the coupling mixes positions and momenta 
variables. 



Remark 3.1 In formula A3.29\) we can see that ifCov x ^(ipo) = then p ur (0) can decrease or increase 
with the coefficient d\ compared to the case d\ = 0. 

Remark 3.2 In the above computations we have assumed that the Planck constant h is one. We can 
rewrite formula ^3.27\ l including h in Weyl quantization and we find 



Pur(0) = -ih 2 tr R 2iv (G T Cov ws GCo V7rE ) + tr K(2d J S GJ E G T ) 



(3.30) 



As it is expected the "decoherence time" becomes very small when h—>0. The quantum correction has 
a meaning only when H is not too small. 
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Remark 3.3 The environment state p T .E is non-negative if and only r < 1. It is enough to prove 
that for N = 1. By an holomorphic extension argument in t we can see that p TjE is negative on the 
subspace of L 2 (R) spanned by the odd Remote functions V>2fe+i where ( — ttt + x 2 )i()k = (2fc + l)ipk- 
Let us prove this. Denote f3 — /3(t) = argtanh(r), and H osc = ^ (— + x 2 ^j . 
By the Mehler formula the Weyl symbol of e~^ Hasc is 

1 -(tanh(ff/2)(x 2 +; 2 ) 

cosh(^/2) 

The thermal state at temperature T — jj is defined as 

■J- _ _^ 

Its Weyl symbol is 

Tpfat) = itanh(/3/2)e- tanh (' 3 / 2 )^ 2 +« 2 ). 
Let us denote Wk(x,^) the Wigner function ofipk- So we have 

(Vfc.e-^-V*) = I e-^ h ^ 2 +^W k (x,OdxaX = c k +^ (3.31) 

cosh(/3) J R2 

Using an holomorphic extension in the variable r, for t > 1, we have 

/3(r) = log J + Y 

Hence we get 

(^Pr,M = ^ 1 1)fc+ ; , Vr > 1. (3.32) 
So if t > 1 jO Ti E lis negative on tpk for k odd. □ 

We can apply these results to give a simple proof that the state the system (S) is corollated with the 
environment (E) for any t > small enough. Here we assume that both /5s (0) and p E (0) are Gaussian 
pure states. A more general result will be given in Proposition 14.151 

Corollary 3.4 Let p(t) = U(t) (p s (0) <S> p E (0))U*(t) be the time evolution of(S)U(E). Assume that 
Purity 7^ (see formula \S.21^ ). Then there exists e > such that for every t e]0,e] we have 

p{t)^Ps{t)®p E (t). (3.33) 

Proof. We shall prove (|3.33p by contradiction. Assume that there exists a sequence of times t n > 0, 
lim t n — such that 

p(tn) = Ps(tn) ® pE(tn)- 

Let us remark that if such a decomposition exists then necessarily we have ps{t n ) — tr E (p(t n )) and 
PE(t n ) = tvs(p(t n )). 
For t — t n we have 

tr(/5(t) 2 ) =ti- s (p S (t) 2 )tT E (p E (t) 2 ) 

From the above result applied to (S) and (£7) we have 

tr(/5(i) 2 ) = (1 - c s i 2 + 0(t 3 ))(l - c E t 2 + 0(t 3 )) = 1 - (cs + c s )i 2 + 0(i 3 ) 

where all the constants cs,c E are positive. So we get a contradiction because tr(/S(t) 2 ) = tr(/5(0) 2 ) is 
independent on t. □ 

In [5] the authors proved for the quantum Brownian motion that is possible to find a Gaussian initial 
state of the system and a temperature T of the environment such that p(t) is not entangled for all 
times in [0, +oo[. 
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4 The master equation in the Weyl-quadratic case 
4.1 General quadratic Hamiltonians 

In this section we find a time dependent partial differential equation (often called the master equation) 
satisfied by the reduced density matrix p$ of the open system (S). Moreover this equation can be 
solved explicitly using the well known characteristics method. 

We extend here to any quadratic Hamiltonians with arbitrary bilinear coupling, several results proved 
in many places |16ill4j[T3lllO) for the quantum Brownian motion model with bilinear position coupling. 
Our results are inspired by the paper [10] but do not use the path integral methods as in the papers 
quoted above. Another difference is that in our case the number N of degrees of freedom for the 
environment is fixed and finite and the number d of degree of freedom for the system is also finite and 
arbitrary. 

There exist many papers in the physicist literature concerning exact or approximated master equation 

for the quantum Brownian motion model (see the Introduction and References in |10|). 

Here we assume that Hs and He are quadratic Weyl symbols of Hs and He respectively. We denote 

$q the Hamiltonian flow generated by Hq and $* the Hamilton flow generated by H = Ha + V. V 

is a bilinear coupling between the system and the environment : V(z, u) = z ■ Gu where G is a linear 

map from M. 2N in R 2d (note that the coupling can mix positions and momenta variables). 

Recall that M. 2d (resp. M. 2N ) is the classical phase space of the system (resp. the environment). 



¥ = $q $* is the interacting flow in the global phase space 



»2(d+7V) 



system is identified to the direct sum 
by 4 matrix blocks. 



The "free" evolution is diagonal 



T>2d , 



p27V 



The phase space of the global 



and in this decomposition the flows are represented 



X]/*. 



It 



$* c 







and the interacting classical Hamiltonian is 



Hi{t, z, u) = V(t, z,u)=z- G(t)u, where G{t) = ($^) T • G ■ 
The classical interacting evolution is given by the equation 



^* = JV 2 u V{t)^\ with J 



Js 
J E 



where V 2 . u V{t) is the Hessian matrix in variables (z, u) e M. 2d ^ n- 
So for the block components of the interacting dynamics we have 



t2N 



¥ u =j s G(t)& ei , ¥ ie 
¥ ei =j E G(t) T %, ¥ e 



J s G(i)** e 
J E G(t) T ¥ ie 



(4.34) 
(4.35) 



Because all the Hamiltonians considered here are quadratic (eventually time dependent) they generate 
well defined quantum dynamics in Hilbert spaces L 2 (R") where n — d (system (S)), n — N (environ- 
ment (E)) and n = d + N (global system (S) U (E)). 



Recall the notations U(t) 



-itH 



, U (t) 



-UHn 



Us(t) 



-itH s 



u E (t) 



-itH. 



E . We have 



Uo(t) = Us(t) ® U E (t) and the quantum interacting dynamics Ui(t) = Uo(t)*U(t). 

At time t = we assume that p(0) = ps(0) ® Pe(0) where ps (resp. p E ) is a density matrix in the 
Hilbert space H s = L 2 (R d ) (resp. in H E = L 2 (R N ). 

p, ps, Pe are the Weyl symbols (i.e the Weyl-Wigner functions, with K = 1, of the corresponding 
density matrices). 
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The coefficients of the quadratic forms Hs, He,V may be time dependent. In this case U(t) means 
U(t, 0) where U(t, s) is the propagator solving 

i^-U(t,s) = HU(t,s), U(s,a)=l H . 

It is well known that for quadratic Hamiltonians, every mixed state p of (S)U(E) propagates according 
the classical evolution (see [6] for details) 

p{t) := U{t)p(0)U(t)*, p(t,X) = p(p,§-*X), X G R 2 ( rf+JV ) (4.36) 

Our aim is to compute ps(t) = trs(/5(t)). The Weyl symbol of ps(t) is given by the following integral 

ps(t,z)= / p(t,z,u)du (4-37) 
Jr 2N 

In particular if p(0) is Gaussian in all the variables X = (z, u) G R 2d x M. 2N then ps(t, z) is Gaussian 
in z. Nevertheless a direct computation on the formula (|4.37p seems not easy to perform, so a different 
strategy will be used. 

4.2 Gaussian mixed states 

Let us consider a very useful class of matrix densities with Gaussian symbols; they are called Gaussian 
mixed states. For the reader convenience we recall here some well known results. 

Definition 4.1 A density matrix p in the Hilbert space L 2 (R") is said Gaussian if its Weyl symbol p 
is a Gaussian p r,m where T is the covariance matrix (positive-definite 2n x 2n matrix) and m G R 2d 
the mean of p. So we have 

p r > m (z) = c r e-i ( - z - m) - r ~ 1{z ~ m) 1 z = (x, G R 2 " (4.38) 

where c 7 = det(r) -1 / 2 . 

Gaussian density matrices are parametrized by their means m and their covariance matrices T where 
we have, as usual 

m= zp T,m (z)dz, ^j,k— / ZjZ k p v,m (z)dz. 

As we have seen in Remark 13.31 some condition is needed on the covariance matrix T such that (I4.38|) 
defines a density matrix (i.e a non negative operator). This condition is a version of the Heisenberg 
uncertainty principle (see a proof in the Appendix). Here J is the matrix of the symplectic for on 
R" x R" and X = (x, £) G R 2 ". 

Proposition 4.2 p T < m (X) defi nes a density matrix if and only ifT + i^ > 0, or equivalently, if and 
only if the symplectic eigenvalues of T are greater than | . 

Moreover, up to a conjugation by a unitary metaplectic transform in L 2 (M. n ), p r '° is a product of one 
degree of freedom thermal states Tp (Remark \3.3\) . 

p ' m {X) determines a pure state if and only if2T is positive and symplectic or equivalently 2T = FF T 
where F is a linear symplectic transformation. 

From these results we can compute the purity (hence the linear entropy) and the von Neumann entropy 
of Gaussian states. 

For the purity we have the straightforward computation: 

PuAp T ' m ) - (2<T n / p F ' m (z) 2 dz = 2-"(detr)- 1 / 2 . (4.39) 

jR 2d 
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Concerning the von Neumann entropy we begin by the computation for thermal states with n = 1 
and temperature T = j . For simplicity we compute with the Neper logarithm "In" . We have 

where r = tanh(/3/2) and ln(e~^ ffosc ) = —(3H osc . Then we have 

-ti(p T lnp T ) - JL-^triHos^ 6 -) + ln(Z(/3)) 
Recall that Z{0) = 2 s inh(j3/2) • ^° we com P u te 



tr(ff osc e-^-) = - Atr(e-^H - -^L 



d/3" y ~ ' 2(sinh/3) 2 
and we get 

S ™M = 2iJb> +b(eW2 - e ~" 2) - 
With the parameter r we get finally the formula: 

8 -^-M^{^)-^)) (4 - 40) 

The parameter r is related with the linear entropy: 

Se(p T ) = l-^l = l-T, (4.41) 

These formulas were obtained in pQ. 

Let us denote £(t) the r.h.s in (|4.40|) . Using Proposition 14.21 and additivity of the von Neumann 
entropy, we get the following formula for a general Gaussian state with covariance matrix T: 

s BN (p T n= E ( 4 - 42 ) 

l<j<n 

where n,-- • ,r„ are the symplectic eigenvalues of (see Appendix for more details). 
4.3 Time evolution of reduced mixed states 

Here we state and prove the main results of this section. It is convenient to work in the interacting 
setting. Recall that we have 

P { s\t) = U* s (t)Ps(Ws(t)), (4.43) 
pf{t,z) = psfr&sz) (4.44) 



Theorem 4.3 Let t c > be the largest time t c such that is invertible for every t G [0,£ c [. Assume 
that the environment density matrix pe is a Gaussian with means tue- Then there exist two time- 
dependent 2d x 2c? matrices A^'(t), B^(t), and a time dependent vector v^\t) G R 2d such that for 
t G [0, t c [ the Weyl symbol p J {t, z) of the interacting evolution pa \t) of the system (S) satisfies the 
following master equation (Fokker- Planck type equation): 

^P { s\t,z) = (A^(t)S7 z )-zp^(t,z) + (B^(t)V z )-Vj s I \t > z) 

+v ( - I) (t)-V zP i p(t,z). (4.45) 
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Moreover we have the following formula to compute A^'it), B^\t) 

A^(t) T = -JsGmt^)- 1 (4-46) 
BW(t) = £ (*) +I,( * )T , where (4.47) 

L(t) = J s -G{t)^l e -^\ l )-H\ c )Cov PE {*\ e ) T (4.48) 
v {I \t) = -JsG(t)& ee m E (4.49) 

Proof. In this proof (and only here) we shall erase the upper index (J) for the interacting dynamics. 
Taking the partial trace in the equation (jl.3p we have 

ps(t,z) = [ {V(t),p(t)}(z,u)du 

For simplicity we shall assume that raj = 0. It is not difficult to take this term into account. 

The Poisson bracket is in the variables (z, u) but due to integration in u we have only to consider the 

Poisson bracket in z. Hence we get 

Ps {t,z) = - [ J s G(t)u-V zP (t,z,u)du (4.50) 
where p(t, z, u) = p(0, ^^(z, u)). 

Denote by /(C) the Fourier transform of / in the variable z. Then we have 

p r s (tX)=i(- [ J s G(t)up{0,^- t (z,u)e- iz< dzdu 

Now let us perform the symplectic change of variable (z',u r ) = ^^ l (z,u). Then we get, using the 
splitting assumption at t = 0, 

hit, = K ■ I dz'du'JsGit)^ 1 + ^ e u')ps{z l )p E {u l )e- K ^ z ' + ^ u '\ (4.51) 

J R 2(d+N) 

Let us denote ip — ^^z' + ^\ e u' and using the equality 

iO^r'Vce-^ - (z' + (^-^tu^e-^ 

we get 

~ps(t,Q = -(■JsG(t)¥ ei (<f t ll )~ 1 VtPs(tX) 

-it- f dz'du'JsGml^r^U^PEiu^psiz^e- 1 ^ 

+<• / dz'du'JsGit^^pE^psiz'^-* (4.52) 

7 R 2(d+N) 

To absorb the linear terms in u' we use that pe is a Gaussian, pe = CAe _1 / 2 "' A ", where A is positive- 
definite 2N x 2N matrix and c\ a normalization constant. So we have A -1 V u p£(it) = upe{u), and 
integrating by parts we have 

Ps{t, = -C • A^(t) T V c p s (t, ()-((■ L(t)()ps(t, 0- (4-53) 

We get (|4.45p by inverse Fourier transform. 
□ 

We can deduce a master equation for the state psif) of the system (S): 
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Corollary 4.4 With the notations of Theorem \4-3\ the reduced density matrix for the system satisfies 
the following master equation 

§- f Ps{t,z) = {H s ,ps(t)}(z) + (A(t)V z )-zps(t,z) + {B(t)V z )-V zPs {t,z) 

+v(t) ■ V zPs {t, z), for t € [0, t c [. (4.54) 

where A(t), B(t) and v(t) are given by formulas \4-4^ with , 5 t = <!>' (replacing the interacting 
dynamics by the complete dynamics and G(t) by G(0) = G). 

Proof. Recall that p(t,z) — p^(t, 3><j*(z)jf|. So the change of variables z i— > $5*2: gives easily the 
result. □ 

Remark 4.5 ps{t,z) is of course well defined for every time t € R &wi i/ie coefficients of the equation 
master ^4-45\ ) may have singularity att = t c as we shall see in examples. We shall see below a physical 
interpretation of t c . 

A consequence is that ps(0) 1— > ps(t) is not always a group of operators. 

Let us give a lower bound for the critical time t c — ini{t > 0, det($*J = 0}. 

Denote 7 = ||G|| (it is a measure of the strength of the interaction) and f(t) — \\Q f s \\ + \\&e\\- 

Proposition 4.6 If sup f(t) < +00 then there exists c > such that t c > — . 

t>o 7 

// there exist C, 5 > such that f(t) < Ce St for every t > then there exists C\ £ R such that 

t c > yl°g(^) + G\. The constants c and C\ are independant of 7. 

Proof. It is enough to work in the interaction representation. Using interacting time evolution ^P* for 
the total classical system (S) U (E), we get 



||¥*-I||<7/ f(s)ds + 1 f(s)\\y s -I\\ds 
Jo Jo 

Denote F(t) — J * f(t)ds. Using the Gronwall Lemma and integrating by parts we get 



II** - I|| < lF(t) + f-F(t) 3 e-< F V (4.55) 
From the inequality (|4.55l) we easily get the Proposition. □ 

Remark 4.7 The coefficients of the equation |^.^<5| ) are related with the first and second moments of 
the reduced density matrix ps(t). Let us denote 

m [ p{t)= [ z j p ( - I \t,z)dz, p i j I k ) (t)= [ z J z k p {I \t 1 z)dz. 

jR2d J K 2d 

From j4-45\ ) we get 

m«(t) = -(A(''(() T m' ; '((), where m^\t) = (m[ T) (t),--- ,m { *}[t) (4.56) 
Computing directly from \4-5Cfy we get 

m { p{t) = - I Zj J s G(t)u ■ S7 zP (t,z,u)dudz = ( J s G(0*L m(/) (0))j- ( 4 -57) 



2d + JV 



So using we 9 e ^ 

m ( i I \t) = * t u m < i I \0). (4.58) 



2 If the Hamiltonians are time dependent, we have to replace $ s by $3 ' 
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From H4-56\ ) and \4-58ty we get again \4-4^ - 

Computations of the second moments, using equation \4-45\ ), gives 

^\t) = 2B^(t) - (A^- T (t)^(t) + »W(t)AW.T(t)). (4.59) 
As above we can also compute directly using ^4-50^ and \4-31$ 

(t) = ^'Homy + ^ eM£ (0)(*L) T . (4.60) 
Using that p^(t) is symmetric we get 

A(/)(t) = 4(*^ (/)(0)( ^ i)T + *^ B(0)(Vl/L)T ) (4 ' 61) 
V {I) (t) = ^(^ (/) (0)(^) T + *L^(0)(*L) T ). (4.62) 

where /i_e(0) is the second moments matrix of pe(0)- 

Remark 4.8 We can get again the formula {3.27ty from \4-45ty - We have ps(0) — and 

p s {0) = A^(0)V zPs ■ z Ps + £«(0)V z p s • V zPs . 

So we get 

i (/) (0) = -J S GJ E G T , flW(0) = -J s GCov PE G T J s 

and 

trs(/Ss(0)/3 s (0)) = itr KM (AW(0)) + (27r)- d / B«(0)V, Ps • V 2 p s dz, 
which coincides with formula \'3.21^ in this particular case. 

Remark 4.9 Suppose that the initial state of the total system (S) U (E) is Gaussian: 

p(0,z,u) = cre"^/ 2 -* 1 " ( Z > U H Z > U ) ; Cr > is a normalization constant. Then we have the following 

direct computation for ps(t, z). 

We get first the Fourier transform: 

w(U)=e- (1/2)rWK ' 0HC,O) 
where T(t) — ($')r($*) T . Using inverse Fourier transform we have 

Ps (t,z) = cr, t e-^S 1 ( t >* (4.63) 

where Ts(t) is the matrix of the positive-definite quadratic form £ i— > r(i)(£, 0) ■ (C>0) on R 2d and 
c r ,t = (2 7 r)- d det 1 / 2 r t . 

We can see on this example what is the meaning of critical times t c . 
Assume that T = T$ © where T$,e are positive-definite quadratic forms on M. 2d respectively M. 2N 
and let us introduce the following quadratic forms on M. 2d . 

C t (C = r((^) T C, (<^) T C) • ((<^) T C, (<^) T C) 

We have Q (C,C) = T S C • C and 

Q t (C, = r^rc • ($*J T C + r E (^) T C ■ (<^) T C (4.64) 

We see from fr4-64\ ! that the initial state ps(0) cannot be recovered from its evolution at time t c : only 
the restriction of Qo to (ker^a) 1 - is recovered from Qt c . The physical interpretation is that a part of 
information contained in ps{tc) has escaped in the environment represented here by Te- 
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We shall see now that the master equation (|4.45[) can be easily solved by the characteristics method 
after a Fourier transform. As in the proof of the Theorem 14. 31 p<j (t, C) denotes the Fourier transform 

Theorem 4.10 Assume that the means of the environment state is 0. With the notations of the 
Theorem \4-3\ we have: 

4 /) (^C) = Ps(o,(^) T C)Gxp(^-ic-e( / )wc) , vteR, (4.65) 

where 0^(t) = *lCov pE(0) (*i) T . 

Remark 4.11 The interpretation of the r.h.s in formula l{4-65\ ) is the following: the first factor is a 
transport term. The second term is a dissipation term due to the influence of the environment which 
is controlled by the non negative matrix Q^'(t). 



Formula {4-.l(ty is an extension of a formula obtained in JlOf for the Brownian quantum motion model 
where Q^\t) is named the "thermal covariance" because the environment is in a thermal equilibrium 
state. 

Proof. In the proof of the Theorem 14.31 we have seen that we have A^(t) T = -*i* ,0) *l°'* ) . 

Here it is convenient to use the notation \In*> s ) for the evolution of ^ at time t starting form initial 

data at time s, as it is usual for non autonomous time dependent evolution equations. 

So the characteristic system for the equation ()4.53|) is the linear differential equation 

C = -(^ 0) v^VC, C(Q) = Co- (4.66) 
It is clear that the solution of (|4~66l) is ((t) = (tf^'^Co- So we have 

p ( P(t,C) =P5(0,(*y T C)exp QfV^Ods) (4-67) 

for < s < t < t c , where b(s, Q = C • 5 (/) (s)C- 

But an easy computation shows that we have the following relation 

B (I \ S ) = \ (A^(t) T e^(t) + e«(t)A«(t) + jO (I) (t) 

Integrating by parts gives 



b ( S ,^ s yods = -±(-6^(tK 



and formula (|4.65|) follows. 
□ 



It is not difficult to go back to the evolution of ps using a change of variable and to the time 
evolution of ps(t) using inverse Fourier transform. 

Corollary 4.12 Under the conditions of Theorem \4-. 10\ we have the following formula 

Ps(t, = ps(0, ($*,) T C) exp (~C ■ 9(t)c) , (4.68) 
with 0(t) = $LCov PE (<i>L) T . 

In particular if p(0) is a Gaussian state with covariance matrix Ts{0) then ps{t) is a Gaussian state 
with the following covariance matrix Tg (t) 

r s (t) = $^r s (o)($y T + $* e r B (o)($L) T , vt e R, (4.69) 

where ^e(0) = Cov PE is the initial covariance of the environment. 
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The formula (|4.68l) is related with Remark 14.71 where we have computed the second moments matrix 
of 

As far as is invertible we see from (I4.68[) that the time evolution of ps is reversible but this is no 
more true for t > t c . 

From formula (|4.68p we get an explicit representation formula for the reduced density of the system 
as a convolution integral: 

Corollary 4.13 With the notations of Corollary \4-12\ we have 

p^(t, z) = c <f>(i) Ps (0, *4°'<V) exp {-\{z - z') ■ ^\t)-\z - z')) dz', 
p s (t,z) = c s (t) ( psiO^^z^expf-liz-z^-Qitr^z-z'^dz'. (4.70) 

where c^(i) = (27r)- d / 2 det(** i )- 1 det (e( ; )) _1/2 and 
c s (t) = (2^)- d / 2 det($* l )- 1 det(6)" 1/2 . 

I/det($* i ) — or if det((<i>* e ($' e ) T ) ~ 0, integrals in \4 ■ ?@\ ) are defined as convolutions of distribu- 
tions. 

Formula (|4.70j) shows clearly the damping influence of the environment on the system because 6(i) 
is a positive matrix under our assumptions. If Q(t) degenerates then the Fourier transform of 
exp (— i£ ■ O(t)C) is a distribution supported in some linear subspace Vs of R 2d where the damp- 
ing takes place. 

We can deduce an exact formula for the linear entropy Si(t) if ps(0) is a Gaussian state. 
First, using the Planchercl formula we have 

S e (t) = l-(27r)- d / (p s (t,z)) 2 dz 

JK 2d 

= l-(2n)- 3d [ \p S (tX)\ 2 d( (4.71) 

JR 2d 

= l-(27r)- d / cxp(-C-e(t)C)|p5(0,(^) T C)| 2 rfC (4.72) 
If ps(0) is a Gaussian state we have 

ft(i) = 1 - det (r s (t)r 1/2 (4.73) 
where Ts(t) is given by the formula (|4.69j) . 

Remark 4.14 We have seen in this section that we can compute the quantum evolution of a reduced 
system when the classical evolution of the blocks $*j and <I>* e of the total system (S) U (E) are known. 
It hopeless to get general explicit formulas for these blocks, but this is possible in the particular case 
of two oscillators (Appendix\^j. For a bath of N oscillators, 1 < N < +oo the problem seems difficult 
(see Appendix\Bty. For continuous distribution of oscillators many results were obtained concerning 
the Quantum Brownian Motion model (see references). 

Let us close this section by giving a simple sufficient condition to get correlations between states of 
the system (S) and the environment (E). If p is a state of the total system (S) U (E) we say that (S) 
and (E) are uncorrelated if p = ps <8> Pe- Note that this decomposition is unique and ps,E = t^E.sP- 
In decoherence theory [8] more difficult notions are also considered: separability and entanglement, p 
is separable if there exists a decomposition p = *^^PjPj <8> pf where J2jPj — ^ Pi — 0, Pj, pf are 

j 

pure states respectively in "Hs, He- If p is not separable it is said that p is entangled. 
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Proposition 4.15 Assume that p(0) = ps(0) ®/3e(0) and that /3s(0), Pe(0) are Gaussian states with 

covariance matrices Ts,Te- 

If the following generic condition is satisfied 

JsTsGJeTe 1 + G £ (4.74) 

then there exists e > such that for every t s]0, e], p(t) ^ Ps(t) <8> f>E(t), which means that the reduced 
system (S) is immediately correlated with its environment (E) when a coupling is switch on. 

Proof. The proof is elementary. 

We compute with the classical Gaussian state in the interaction representation (as above for simplicity 
we erase the index (J)). We have 

p(t, z, u) = p s (0, ^z + *L°' t) U )p B (0, + 9%%). 

Using the assumptions, we have 

p(t, z, u) = K(t, z)L(t, u)e- z - c ^ u (4.75) 

where the crossed term is determined by the matrix C(t) satisfying C(0) = 0. Using the interacting 
evolution equation we can easily compute the derivative (7(0). So we find C(0) = T^ 1 J$G — GJe^e 1 
and C(0) ^ if condition (14.741) is satisfied. Hence C(t) ^ for t small enough and t ^ □ 

Remark 4.16 In J8$ the authors have proven a stronger entanglement result for the Quantum Brow- 
nian Motion model. 



A Coupled Harmonic oscillators 



More explicit computations can be done for systems with two coupled one dimensional harmonic os- 
cillators. 

Hs = \{-d 2 x , + ^ s V), H E = \{-0 2 y2 + uo E 2 y 2 )- 



We assume that uig > and uje can be real positive or purely imaginary (lu e < 0). lo\ > means a 
stable environment and uj 2 e < means an unstable environment. Unstable environment was consid 
ered in the paper 0]. Here we revisit the computations of [I]. 



Le us compute the classical flow of the total Hamiltonian 

H(x, £, y, n) = i(£ 2 + if + u s 2 x 2 + oj E 2 y 2 ) + jxy. 

Let M = ^ 2 ^ be the matrix of the quadratic potential for the total system: V(x, y) 

ojs 2 x 2 + uiE 2 y 2 + jxy. 

The eigenvalues A± of M are 



So we have 



A l = 2 ( ~'e '< yjui-^l)- + 4-- 
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where P = 



cos 9 
sin 9 



— sin 9 
cos 9 



and 



cos 9 



sintf = 



4 



1 Ug - Ul 

2 ( + (c4-c4) 2 + 4 7 2 ) 1 /2 



(1 



1/2 



1/2 



(A.76) 



Now we compute the matrix of $' in the canonical basis of M. x x 

/°\ 

_ 1 

' 6y ~ 
V o / V o J 

IT -1 where 



e x = 




1 

v o y 






V i J 



We have $* = n$^ iag 



P 
2 



o 2 
p 



and 



$diag 



cos(^Ay 2 ) 


-X 1 / 2 sin(tA+ /2 ) 




n = 
o 

cos(tA 


-Ai /2 sin(tAi /2 ) 



'1/2 



l/2x 



sin(tA+ /2 ) 




-1/2 



cos(tA+ 2 ) 






sin(tA 


cos (a V2 ) 



l/2x 



The coefficients of the matrix $* are easily computed using the formula 

e a ■ = n _1 e Q • s^n-^. 

Let us denote k) the coefficient with raw index j and column index k of the matrix $*. This 
matrix has four 2x2 blocks: ^\ i ,<^\ e ,<^ t ei ^ t ee . (16 terms). So we have <&*j(j, fc) = $ l (j,k) for 

e {1,2}, $f e (j,fc) - with j e {1,2}, fc e {3,4}, $* ei (j,fc) = &(j,k) with j e {3,4}, 

fc G {1, 2}, $* e (j, k) = fc) with j, k G {3, 4}, 



$*,(1,1) = cos^cos^A^ + sin^cos^A 172 ) 

$* 4 (1,2) = A; 1/2 cos 2 0sin(a+ /2 ) + A: 1/2 sin 2 0sin(tA 1/2 ) 

$*i(2,l) = -Ay 2 cos 2 0sin(a+ /2 ) - A^sin^sin^A^ 2 ) 

$*,(2,2) = cos 2 61 cos (i A 1 / 2 ) + sin 2 6 1 cos^A 172 ) 



We have analogous formula for the blocks $- e , & ei , $* e . 



$1(1,1) 
$1(1,2) 

*ie(2,l) 

$1(2,2) 

We note that $* ei = 

$L(M) 

$L(1,2) 
$L(2,1) 
$L(2,2) 



sinflcos^cos^A^ 2 ) — cos(iA_)) 



sin 9 cos 9(X + 1 ^ 2 sin(i A 



,-1/2 



sin^A 172 )) 



sin 6» cos 9(-X 1 / 2 sin(iAi /2 ) + A 1/2 sin(f A V2 )) 



sin6'cos6 1 (cos(iA 1 ( 72 ) — cos(fAi // )) 

cos 2 9 cos(a V2 ) + sin 2 6» cos(iA+ /2 ) 
AI 1/2 cos 2 9 sinitX 1 / 2 ) + A; 1/2 sin 2 9 sin(af ) 
-A 172 cos 2 6> sin(a V2 ) - X 1 / 2 sin 2 6> sin(tA+ /2 ) 
cos 2 9 cos(a V2 ) + sin 2 6» cos(tA+ /2 ) 



vV2^ 



(A.77) 
(A.78) 
(A.79) 
(A.80) 

(A.81) 
(A.82) 
(A.83) 
(A.84) 

(A.85) 
(A.86) 
(A.87) 
(A.88) 
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In particular we have 

det =1 + 2 sin 2 6 cos 2 8 
(cositX 1 / 2 ) cos(U V2 ) + X ~ sin(tA+ /2 ) sin(ai /2 ) - lj (A. 



and 



det(^)=l + 2- 



( W 2 ~^) 2 + 4 7 2 



( cos(tA + /2 ) cos(U 1/2 ) + fs + <»l ^/Jj sin(tA V 2) _ A (A go) 

From equation (|A.90I) we can compute the critical time t c . We see here that we may have t c finite or 
infinite, depending on the coupling constant 7 and on the sign of u> E . 

1. Assume that oJs^e > 7 2 an d w l uj 'e- 

We have A+ > A_ > 0. We see easily that, for somme constant C depending only on ujs,^e, we have 

I det -1| < C7, VtR 
Then if 7 < ^, ($^) is invertible for every tel. 

2. If (Jg = u>% and w| > 7 2 we get from (|A.90[) that there exists C > such that is invertible for 

< t < e . 
— — 7 

3. If < 7 2 then A + > > A_. We can find constants C\, C 2 such that if 7e C2 * < C\ then $* 4 is 
invertible. 



B On the quantum Brownian motion model 

The quantum brownian motion model was introduced in the paper |12) . The starting point is a 
one degree of freedom system interacting with a large number N of harmonic oscillators. The total 
(classical) Hamiltonian is the following 

H(x,t;y,r,) = £-+V(x)+ ]T ^ + k 3 {y 3 ~x) 2 . (B.91) 

i<j<N J 

where x, £ are the coordinates of the system, y = (y± , • • • , j/jv) an d f] — (771 , • • • , tjn) are the coordinates 
for the environment. The environment consists in N harmonic oscillators and the system is connected 

to each oscillator by a spring with constant kj > 0. This model is a particular case of quadratic systems 

2 

considered in section |H if V(x) = "^§-x 2 . The Hamiltonian can be splitted as H — Hs + He + Hi 
where 



H s = 


e 

2m 


H E = 


»( 


Hi = 


-2x 



V(x) + \{ k >> 2 (B.92) 



^ l<j<N 



E £ • Ml ] (B.93) 

(B.94) 



m, 

l<j<N 3 



l<j<N 



The classical evolution for the Hamiltonian is not explicitly given for N > 1 (we have got explicit 
formula in Appendix lAl if N = 1). It can be seen that the time evolution of the position x of the 
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system satisfies the following equation assuming m 3 - = 1 for simplicity: 

mx(t)+ [ K(t-s)x(s)ds + V'(x(t))+K(t)x(0) = F(t) (B.95) 
Jo 

where 

K(t) = kj cosiVkjt) (B.96) 

l<j<JV 

F(t) = Y Vj (0)cos(Vk jt) + r) j (0)sm(Vk j t) (B.97) 

l<j<JV 

The difficulty here is that equation (|B.95|) is not an ODE because of the integral term. 
In [12l [11] the authors considered a large N limit, and a continuous distribution of oscillators, such 
that (|B.95[) is transformed into a stochastic differential equation where the integral term is replaced 
by the damping term jx(t) where 7 is a damping constant (Langevin equation). 



C Gaussian density matrices 

We shall give here a proof of Proposition ^. 21 This is consequence of a particular case of the following 
Williamson theorem (see [35] and [5] Appendix 6). 

Theorem C.l Let V be a positive non degenerate linear transformation in M. 2N . Then there exists a 
linear symplectic transformation S and positive real numbers X\, ■ ■ ■ , Xn such that 

S T TSej = Xjej, and S T TSej + N — Xjej + N (C.98) 
for 1 < j < N, where {e±, • • • , ejv, ■ ■ ■ , ^2n} is the canonical basis ofM. 2N . 
in [5D] the authors gave a simple proof that we recall here. Recall the following known lemma 

Lemma C.2 Let A be a non degenerate antisymmetric linear mapping in M. 2N . Then there exists an 
orthonormal basis {v%, u*, • • • , ujv, v^} ofM. 2N and positive real numbers {v\, ■ ■ ■ , vpf}, such that 



Avj = VjVj, Av* = 



Proof of Lemma ] C. Si We proceed by induction on N. This is is obvious for N = 1. 

Assume N > 2. Let v\ be an eigenvalue of the symmetric matrix A 2 , A 2 v\ = V\V\, \v\\ — 1. We 

can choose a vector v\ and a € K such that Av\ = —av\ and \\v\\\ = 1. Then we have easily that 

v\ ■ v\ = and if P is the plane spanned by {v%, u*} then P and P 1 - are invariant by A. So we can 

apply the induction assumption to A acting in P and the Lemma is proved. □ 

Proof of Theorem \C.l\ Consider the antisymmetric matrix A = r _1 / 2 jr -1 / 2 . Using Lemma [C.2I we 

can find an orthogonal matrix R and a diagonal matrix — diagjVi, 1/2, ■ ■ ■ , vn}, Vj > 0, 1 < j < iV 

such that 

/q-1/2 n \ 

Denote B=l 1/2 J and S = Y~ 1 / 2 RB. We get easily that S T TS = B 2 and S T JS = J so 



the proof of the Theorem I C . 1 1 follows . □ 

The N real numbers Xj in Theorem I C . 1 1 are the symplectic eigenvalues of T. Using that S T = —JS~ 1 J 
we see that JT is diagonalizable with eigenvalues {±Aj, 1 < j < N} and that {Ai, • • • , Ajv} are the 
eigenvalues (with multiplicities) of \ JT\ = (— JT 2 J) 1 / 2 . 
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Proof of Proposition ^. 2\ We can assume assume that m = and we denote pr = p r ' - 

We use the symplectic normal form for F given by Theorem lC.il Let R(S) be the metaplectic unitary 

operator associated with 5" (see for example [5]). Hence we have 

R(S)*p r R(S) = p s rrs = Pn ® p T1 ■ • • ® Pt n , 

where the Tj are the symplectic eigenvalues of ■ So pr is a density matrix if and only if we have 
< Tj < 1 for 1 < j < d (Remark 13. 3|) . This condition means that the symplectic eigenvalues of 2T 
are greater than 1 or equivalently that 2T + i J > 0. 

We have already seen that if pr is a pure state then we have 2T = F T F with F symplectic. Conversely 
if 2r = F T F then pr is the Wigner function of a squeezed state R(F)ipo, (po being the standard 
Gaussian (for details see [6]). 

Gaussian states are thermal states for positive non degenerate quadratic Hamiltonians and con- 
versely. This can be proved as follows. 

Let H (it) = \u ■ Mu where A is a positive non degenerate linear transformation in R 2Ar . We consider 
symplectic coordinates for it: yj — u ■ ej and rjj = u ■ ej+N, 1 < j < N. 

Using Theorem IC.ll we have H(u) = Ha{Su) where S is a symplectic linear transformation and 
Ha{u) = Aj- u.^ , Uj = {xj,£j) £ M 2 . Applying the metaplectic transformation we have H = 

l<j<N 

R(S)H A R(S)* and 

e- pk =R(S)e-^R(S)* 
From the Melher formula for the harmonic oscillator the Weyl symbol Wa of e~^ H& is 

tit ( \ TT cxp(-tanh(/3A 3 )(y 2 +ry 2 )) 

WaM= M» — — 

so the Weyl symbol W of e~@ H is given by 

W(y,r,) = W A (S- 1 (y,v)). 
This proves that — l> _ f}iI is a Gaussian state. More explicit results are given in [71 115). 
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